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ABSTRACT 

The synergy between the Fermi-LAT and ground-based Cherenkov telescope arrays gives us the 
opportunity for the first time to characterize the high-energy emission from blazars over 5 decades in 
energy, from 100 MeV to lOTeV. In this study, we perform a Fermi-LAT spectral analysis for TeV- 
detected blazars and combine it with archival TeV data. We examine the observational properties in 
the 7-ray band of our sample of TeV-detected blazars and compare the results with X-ray and GeV- 
selected populations. The spectral energy distributions (SEDs) that result from combining Fermi-LAT 
and ground-based spectra are studied in detail. Simple parameterizations such as a power-law function 
do not always reproduce the high-energy SEDs, where spectral features that could indicate intrinsic 
absorption are observed. 

Subject headings: galaxies: active - galaxies: nuclei - gamma rays: general 
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1. INTRODUCTION 

Active Galactic Nuclei (AGNs) are extreme objects 
with observed luminosity outshining their host galaxy. 
These sources are believed to be powered by accretion 
onto a central supermassive black hole, commonly dis- 
play relativistic jets, and exhibit non-thermal contin- 
uum emission extending from the radio band to X and 
7 rays. Blazars constitute a subclass of AGNs, with jet 
axes oriented close to the observer's line of sight. Rel- 
ativistic beaming gives rise to distinctive observational 
features in blazars, such as strongly anisotropic radia- 
tion, super luminal motion, high p olarization and rapid 
variability (jUrrv fc Padovamlll995D . Blazars are divided 
into two subclasses, flat spectrum radio quasars (FS- 
RQs) and BL Lacertae objects (BL Lacs). FSRQs are 
observationally characterized by broad spectral lines in 
the optical band, which are weak or not present in BL 
Lacs. The spectral energy distribution (SED) of blazars 
exhibits a two-component structure, with a low-energy 
component peaking between infrared (IR) and X-ray en- 
ergies, and a high-energy one between X and 7 rays. The 
low energy component is believed to be dominated by 
syn chrotron emission from rela tivistic electrons in the 
jet (|Kembhavi fc Narlikarl [199 9). The peak frequency of 
the synchrotron component of the SED (24 y n) is used to 
sub-classify BL Lacs into low (LBLs, v syn < 10 14 Hz), 



intermediate (IBL, v syn 
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10 15 Hz) and high- 



frequency-peaked BL Lacs (HBL, v syn > 10 Hz). 

The high-energy component of the blazar SED has 
been historically less studied, due to the later devel- 
opment of hard X-ray and 7-ray detectors compared to 
those of longer frequency bands. The Synchrotron self- 
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Compton (SSC) model is the simplest scenario that ex- 
plains the high-energy emission of blazars, by inverse- 
Compton (IC) up-scattering of soft synchrotron photons 
off the s ame electrons that ha ve undergone synchrotron 
cooling (iMaraschi et al.1 fl992l ). Throughout the text, 
we refer to the high-energy component of the blazar 
SED as "IC component" . An additional IC target pho- 
ton field external to the jet is often invoked (see, e.g., 
IW Comae: Acciari et al . 2009d). This mechanism is re- 
ferred to as external-Compton (EC), and several pos- 
sible sources for the external photon field have been set 
forth (|Bottcheiil2010f ). Other models suggest a significant 
contrib ution from hadronic processes to the high-energy 
output (Man nheim fc B iermami 1993). 

A good spectral characterization of the high-energy 
peak of the blazar spectrum (keV - TeV band) 
is essential to discriminate between the aforemen- 
tioned mo dels. During the EG RET era covering 
1991-2001 (jThompson et al.l Il993h . only five blazars 
were known at TeV energies: Mrk 421, Mrk 501, 
1ES 1959+650, PKS 2155-304 and 1ES 2344+514; 
thanks to the first generation of ground-based Imag- 
ing Atmospheric Cher enkov Tele scopes (IACTs; Whip - 
ple (jKildea et al.ll2007h. HEGRA (IPuhlhofer et aUjflOg) , 
Durha m Mark 6 ([Armstrong et al.lll999l ). and Telescope 
Array (|Aiso et al.lll997f) V Only three of these sources 
were detected by EGRET (Mrk 501 was only marginally 
detected and 1ES 2344+514 not seen at all). By the 
time Fermi started operations (2008 August), the num- 
ber of known TeV blazars had increased to 21 with 
the s econd generation o f IACTs i n operation (VERI - 
TAS (IWeekes et alJl20Toh. M AGIC (lAleksic et al.ll20rl . 
HESS (|Bernl5hr et alj l2003)). This number has doubled 
since then, with most TeV blazars being also detected in 
the GeV range by Fermi-LAT . 

For the first time now, good quality spectra are avail- 
able both from Fermi-LAT in the high-energy (HE, 
0.1GeV< E < 100 GeV) 7-ray band and IACTs in the 
very high energy (VHE, E > 100 GeV) 7-ray band for 
more than two dozen sources. The combined spectral 
data covers up to five decades in energy giving a detailed 
description of the high-energy peak of the blazar SED. 
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Fig. 1. — Skymap of TeV blazars in galactic coordinates, as of 2012 January, generated using TeVCat (http://tevcat.uchicago.edu/ ). 
Blue and pink shaded areas represent VERITAS/MAGIC and HESS visibilities, respectively. A total of 46 sources consisting of 33 HBLs, 4 
IBLs, 4 LBLs, 3 FSRQs, and 2 sources that were formerly c lassified as AGN of un known type (UNID), namely IC 310 and VER J0521+211, 
are shown. VER J0521+211 is now identified as a BL Lac (Err ando et al.H2011bl) a nd recent studies suggest that the high-energy radiation 
from IC 310 originates from a blazar-like emission mechanism IIKadler et al.N2012i) . 



Recent studies have explored this newly available data 
sample, focusing on the G eV properties of TeV-selected 
blazars (|Abdo et al.l I2009D , or deri ving jet parameter s 
assuming leptonic emission models (|Zhang et all 120121) . 
These studies are similar to earlier studies carrie d out on 
a lim ited sample of TeV-detected blazars fe.g.. IWagnerl 
12001 , 

In this paper we study the GeV- TeV observational 
properties of the high-energy emission in blazars that 
are detected in the TeV band. Section 2 describes the 
population of TeV blazars, giving census information, 
investigating luminosity, redshift and photon index dis- 
tributions among different blazar types. In Section 3, 
we study TeV blazars that appear in the Fermi data 
and outline their GeV properties with respect to the rest 
of the Fermi blazars. Section 4 defines our sample and 
focuses on general observational TeV properties of our 
objects. In section 5, we give a detailed description of 
the Fermi analysis that we performed on our TeV blazar 
sample. Finally, Section 6 discusses various observational 
characteristics of the studied sources based on their GeV- 
TeV spectral shapes, such as the peak frequency of the 
IC component, absorption-like spectral features and vari- 
ability. Throughout the text, the symbol a is used to des- 
ignate the standard deviation, as a measure of statistical 
significance. 

2. TEV BLAZARS 

Mkn 421 was the first blazar and extragalactic object 
to be discovered as a VHE 7-r ay emitter, detect ed with 
the Whipple telescope in 1992 ([Punch, et al.lll992t ). Since 
then, different candidate selection methods have been ap- 
plied to radio, X-ray or HE data with the aim of finding 
new "TeV" blazars, i.e. detected in the VHE regime, 
dCostamante fc Ghisellinil 120021 : Ide la Calle Perez et al.l 
120031 : iBehera fc Wagnerl 120091 ). leading to the discovery 
of most of the known TeV blazars. To date, 44 blazars 
and 2 AGNs of unknown type have been detected in the 
VHE range0, with a census consisting of 33 HBLs, 4 IBLs, 

6 http://tevcat.uchicago.edu/ 



4 LBLs and 3 FSRQs (see Figure [T]). In this work, we 
have studied the blazars that have a published TeV spec- 
trum as of 2011 February (referred to as the "sample" in 
the remainder of the text). 

The redshift (z) of TeV blazars in our sample ranges 
from 0.031 (HBL Mrk421) to 0.536 (FSRQ 3C279), and 
nearly one fourth of the population does not have a se- 
cure redshift. This lack is due to the fact that optical 
emission lines are typically weak or absent in BL Lac 
objects, rendering direct redshift measurements difficult. 
The majority of known-redshift TeV blazars are located 
at z < 0.2, mostly due to the absorption from the ex- 
tragalactic background light (EBL). Figure [5] illustrates 
the redshift distribution of all TeV blazars. The TeV FS- 
RQs are the most distant objects in the population with 
redshift ranging from 0.36 to 0.536. 

Using archival data, we calculated the apparent 
isotropic luminosity of the blazars in our sample (see Sec- 
tion 0} for E > 400 GeV, with the following formula: 



L = F x AttD 2 l /(1 + z) 



2-r 



(1) 



where F is the energy flux for energies above E > 400 
GeV, Dl is the luminosity distance with Hubble con- 
stant Fl = 71kms~ 1 Mpc~ 1 and the cosmological con- 
stant £l\ = 0.730, z is the redshift, and T is the observed 
photon index for each blazar. Figure [3] top panel shows 
the luminosity versus redshift correlation for the sample. 
Sources at high redshifts tend to be scarce and much 
more luminous. The reason why we see only luminous 
sources at high redshifts is that the less luminous ones 
are too weak to be detected. On the other hand, the 
reason why we do not detect them in low redshifts is 
that we integrate over a much smaller volume and thus 
are less likely to see high-luminosity sources that should 
be scarce compared to low-luminosity ones. Note that if 
the luminosity is corrected for EBL absorption, which is 
stronger at TeV energies and high redshifts, the correla- 
tion will be steeper. Figure [3] bottom panel shows the 
TeV photon index versus luminosity correlation of the 
same sample. For luminosities up to ~ 10 45 ergs _1 , the 
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Fig. 2.— Redshift (z) distribution of all TeV blazars as of 2012 
January. Top, middle and bottom panels also show HBLs, IBLs 
and LBLs, and FSRQs, respectively. Most of the blazars with 
known redshifts (30 out of 39) are located at z < 0.2. FSRQs are 
the most distant objects in the population. The farthest object 
is the FSRQ 3C 279, with a redshift of 0.536. The rest of the 
population does not have a secure redshift. 

photon index distribution is fairly homogenous. 

3. Fermi TEV BLAZARS 

The second Fermi-LAT catalog (2FGL) contains 1873 
sources, among which 1062 are AGN, with 435 BL Lacs, 
370 F SRQs and 257 AGNs of unknown class (jNolan et al.l 
12012ft . Thirty six of these AGNs are TeV emitters. Fig- 
ure |4] (top) shows the distribution of Fermi spectral in- 
dices (rQev) for TeV and non-TeV blazars in the 2FGL 
catalog. TeV-detected blazars tend to have harder GeV 
indices. As can be seen from Figure U (bottom), an- 
other distinguishing parameter for TeV emitters within 
the Fermi blazar population is the integral flux for ener- 
gies above 1 GeV. It follows that an effective method for 
TeV-candidate selection in the HE 7-ray band is to look 
for bright hard spectrum sources, and select the candi- 
dates based on their extrapolated fluxes at VHE ener- 
gies. For all TeV blazars, TgcV < 2.3 and for most of 
them TgoV < 2, in agreement with an inverse-Compton 
peak frequency (j^ic) located in the high-energy tail of 
the Fermi range or beyond. Figure [5] shows scatter plots 
of spe ctral indices of fermi-bright AGN (jAbdo et al.l 
l2010a|) in radio, optical and X-ray bands, comparing TeV 
and non-TeV sources. The TeV and non-TeV AGNs oc- 
cupy separate regio ns in the paramete r space, consistent 
with the results in ( Abdo et al.ll2010aft . considering that 
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Fig. 3. — Top: apparent luminosity (not EBL-corrccted) vs. red- 
shift plot of TeV blazars. As expected, high-redshift objects are 
less likely to be seen, and in fact most of them were detected in flare 
states. Bottom: TeV index vs. apparent luminosity distribution of 
the same sample. For fluxes up to ~ 10 45 ergs — 1 , the TeV index 
distribution is fairly homogenous. At higher luminosities, mostly 
soft spectrum sources are detected. FSRQs are the most luminous 
TeV emitters but one should note that their data come from flares 
only. Squares represent HBLs, triangles IBLs and LBLs, and cir- 
cles FSRQs. Points representing different flux states of the same 
source are connected with dashed black lines (bottom). 

most TeV AGNs are HBLs. 

4. DATA SAMPLE 

Our blazar sample contains all blazars with a pub- 
lished VHE spectrum before 2011 February, including 
a total of 26 sources (see Table [J): 19 HBLs, 3 IBLs, 2 
LBLs and 2 FSRQs. TeV spectral index distributions 
of the whole sample are shown in Figure [6] Three of 
these blazars have insecure redshifts cither because 
the spectroscopic measurements were inconclusive, or 
the calculations were made indirectly based on EBL 
absorption studies. References for the adopted redshift 
values in these three cases are given in Table [TJ Seven 
of our targets were detected with EGRET and 23 of 
them are in the 2FGL catalog (jNolan et al.ll2012D . The 
ones that are missing in the Fermi data (1ES 0229+200, 
1ES 0347-121, PKS 0548-322) are very hard spectrum 
sources that would be weak in the Fermi band. More 
than half of the sample have been detected multiple 
times in the VHE band. These multiple detections 
extending over several years and obtained mostly with 
different instruments suggest that spectral variability 
in the VHE band is a common property for VHE 
blazars. Even though no general pattern has been 
established for VHE variability, several sources have 
been observed to have a flux i ncrease up to a few times 
their baseline emission (|W Comae: Acciari et al.l 
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Fig. 4. — Top: GeV photon index (rQ e v) distribution for non- 
TeV blazars (black solid li ne) and TeV bl azars (red shaded area) 
from the 2FGL catalog (Nol an et alj|2012ft . TeV emitters tend to 
have hard spectra in the GeV band. Bottom: Tq c v vs. integral 
flux for E > lGeV for the same sample, with gray open squares 
for non-TeV BL Lacs, black open circles for non-TeV FSRQs, red 
triangles for TeV BL Lacs and dark red inverse triangles for TeV 
FSRQs. TeV emitters are mostly located along the bottom of the 
distribution. TeV FSRQs are at the edge of the group where the 
luminosity is relatively higher and the index softer. 
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iMkn 501: Albert et al.l 



PKS 2155-304: Aharonian et al 



[2007el 

I2007dl) . occasion- 
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ally accompanied by a 
dex (|Mkn 501: Albert et all . 

flux doubling times (IMkn 501: Albert et all l2007et 
IPKS 2155-304: Aharonian et alJl2007d[ ). 

The first 27 month Fermi data and archival VHE spec- 
tra published before 2011 February were used to con- 
struct combined GeV- TeV SEDs in this study. Only 
in seven cases (RGB J0710+591, 1ES 1218+304, PKS 
1222+21 (4C +21.35), PKS 1424+240, PKS 2155-304, 
and two different measurements for 3C 66A) were the 
VHE data found to overlap with the Fermi era. The re- 
mainder of the VHE data were taken before the Fermi 
mission. 

All VHE spectra were corrected for the EBL absorp- 
tion using the model by (|Dominguez et ai]|2011l). Other 
backg round models are also available (e.g.. iFinke et al.l 
l2010f f. However, with a different EBL model, we do not 
expect any significant differences in our results up to a 
few TeV, given the redshift and energy range of our sam- 
ple. See Section |6~21 for a more detailed discussion on the 
EBL correction effects on our study. 

5. FERMI ANALYSIS 

The fact that most of the GeV and TeV data are not 
contemporaneous makes it hard to interpret the com- 
bined spectra of blazars. Moreover, Fermi data represent 



an average state over relatively long periods, whereas 
the VHE spectra consist of "snapshots", mostly taken 
during flares. To account for blazar variability and 
the non-contemporaneous nature of the data set, for 
bright enough sources, the Fermi data were split into 
"low" and "high" flux states as described below. Thus, 
non-contemporaneous GeV and TeV measurements were 
matched in a more realistic way than directly using all 
the time-averaged Fermi data. Table [1] summarizes the 
Fermi flux states and VHE spectra used for each source. 

For VHE data that were taken during the Fermi era, 
time periods of a few months that cover the correspond- 
ing VHE observations were selected for the Fermi spec- 
tral analysis. For blazars that have VHE spectra mea- 
sured before the Fermi era, the first 27-month of Fermi 
data were analyzed (from 2008 August 4 to 2010 Novem- 
ber 4). In all the analysis steps, an energy selection from 
300 McV to 100 GeV was applied to the data. 

The Fermi data were analyzed in the following way. 
First, a 27-month light curve analysis was performed 
for each source using an aperture photometry technique. 
Diffuse class events from a region of 1° radius from the 
target location were selected and counts were plotted as 
a function of time, each time bin containing 49 counts, 
corresponding to a signal to noise ratio of 7. For sources 
with high statistics, low- and high- flux states were iden- 
tified and separated using the average count rate as a 
threshold. Figure [7] shows the resulting light curves for 
all sources, with fluxes normalized to arbitrary units. It 
should be noted that in this analysis, no background sub- 
traction was performed and therefore the resulting light 
curves merely give an estimate of high- and low-state 
time slices. 

Next, a spectral analysis was done for each data 
set. Diffuse class events from a region of interest of 
8° radius were selected and analyzed with Fermi Sci- 
ence Tools v9rl8pd using instrument response functions 
P3-V6-DIFFUSE Sources from th e first Fermi-LAT 
(1FGL) catalog (|Abdo et alJl2010b[ ). bright spots with 
test statistics > 25 and standard galactic and isotropic 
diffuse emission background component^ within the 
region of interest were included in the source model 
files. Un binned max i mum-likelihood an alysis as de- 
scribed in ||Cashlll979t iMattox et all 119961 ) was applied 
to each data set, assuming a power-law (PL) spectrum 
as given in Equation 



dN/dE = N (E/E y 



(2) 



Additionally, to look for possible spectral features in the 
data, spectral points were calculated and fitted with dif- 
ferent power-law functions, and the results were com- 
pared. See Section RT41 for more details. 

Finally, combined GeV- TeV SED data sets were con- 
structed using archival TeV spectra and the correspond- 
ing flux state information from references shown in Ta- 
bled] With each TeV spectrum, the most suitable Fermi 
data subset (average, low- or high- state) was used for 
further study. 

6. RESULTS AND DISCUSSION 
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Fig. 5.— Scatter plots of o ro (5 GHz - 5000 A) n ox (5000 A- Ike V) and a x7 (IkeV - 100 MeV) spectral indices for the Fermi-LAT 
bright AGN sample (LBAS). Data are taken from (Abdo ct al. 2010a). Note the correlation between a TO and a x7 . Red triangles represent 
TeV-detected AGNs, and asterisks non-TeV-detected ones. TeV-emitters seem to be well isolated in a IO - 0x7 parameter space. 



Name 


SED type 


z 


log 10 (* ayn /lHz) 


Fermi var. 


Fermi state 


r TeV 


I 


E th (GeV) 


Reference 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V 


(8) 


(9) 


RGB ,10152 + 017 


HBL 


0.080 






average 


2.95 ± 0.36 stat ± 0-20 syst 


2% 


300 


(Aharonian et al. 2008) 


3C 66A* 


IBL 


0.444 a 


15.63 


171 


MAGIC 


3.64 ± 0.39 stat ± 0-25 syst 


6% 


200 


(Alcksic ct al. 2011a) 












VERITAS 


4.1 ± 0.4 stat ± 0.6 svst 


8% 


100 


(Acciari et al. 2009b) 


1ES 0229 + 200 


HBL 


0.140 


19.45 




average 


2.50 ± 0.19 stat ± 0.10 svst 


2% 


580 


(Aharonian et al. 2007b) 


1ES 0347-121 


HBL 


0.188 


17.94 




average 


3.10 ± 0.23 stat ± 0.10 syst 


2% 


250 


(Aharonian ct al. 2007a) 


PKS 0548-322 


HBL 


0.069 


16.84 




average 


2.86 ± 0.34 stat ± 0.10 syst 


1% 


200 


(Aharonian et al. 2010) 


RGB .10710 + 591 


HBL 


0.125 


21.05 


6 


VERITAS 


2.69 ± 0.26 stat ± 0.20 syst 


3% 


300 


(Acciari ct al. 2010c) 


S5 0716 + 714 


LBL 


0.300 


14.46 


266 


high 


3.45 ± 0.54 sta t ± 0.2 svst 


9% 


400 


fAndcrhub ct al. 2009a.) 


1ES 0806 + 524 


HBL 


0.138 


16.56 


20 




3.6 ± 1.0 st at ± 0.3 syst 


2% 


300 


(Acciari et al. 2009a) 


1ES 1011 + 496 


HBL 


0.212 


16.74 


16 


high 


4.0 ± 0.5 s tat ± 0-2 sys t 


6% 


200 


(Albert et al. 2007d) 


1ES 1101-232 


HBL 


0.186 


16.88 


1 


average 


2.94 ± 0.20 sta t 


3% 


225 


(Aharonian ct al. 2007c) 


Markarian 421* 


HBL 


0.031 


18.49 


44 


medium 


2.20 ± 0.08 stat ± 0.2 svst 


50-200% 


200 


(Albert et al. 2007b) 


Markarian 180 


HBL 


0.046 


18.61 


10 


average 


3.3 ± 0.7 stat ± 0.2 svst 


11% 


200 


(Albert ct al. 2006) 


1ES 1218+304* 


HBL 


0.182 


19.14 


15 


average 


3.08 ± 0.34 Btat ± 0.2 syst 


7% 


200 


^cciari ct al. 2009c) 












VERITAS 


3.07 ± 0.09 stat 


6% 


200 


(Acciari et al. 2010b) 


W Comae* 


IBL 


0.102 


14.84 


47 


high 


3.81 ± 0.35 stat ± 0.34 syst 


9% 


200 


(Acciari et al. 2008) 


PKS 1222 + 21 


FSRQ 


0.432 


13.27 


101 


MAGIC 


3.75 ± 0.27 stat ± 0.2 syst 


100% 


100 


(Alckaic ct al. 2011b) 


3C 279 


FSRQ 


0.536 


12.67 


898 


high 


4.1 ± 0.7 Bt at ± 0-2 syst 


15% 


200 


(Albert et al. 2008a) 


PKS 1424+240 


IBL 


0.260 b 


15.7 


26 


VERITAS 


3.80 ± 0.5 stat ± 0.3 syst 


3% 


140 


(Acciari et al. 2010a) 


H 1426+428 


HBL 


0.129 


18.55 


7 


average 




3% 


1000 


(Horns ct al. 2004) 


PG 1553 + 113 


HBL 


0.4 C 


16.49 


44 


high 


3.4 ± 0.1 s tat ± 0-2 Bys t 


8% 


200 


(Alcksic et al. 2010) 


Markarian 501* 


HBL 


0.034 


16.84 


46 


low 


2.79 ± 0.12 stat 


20% 


200 


(Andcrhub ct al. 2009b) 


1ES 1959+650* 


HBL 


0.048 


18.03 


16 


low 


2.58 ± 0.18 stat 


10% 


200 


(Tagliaferri et al. 2008b) 


PKS 2005-489* 


HBL 


0.071 




9 


average 


3.20 ± 0.16 stat ± 0.10 syst 


3% 


400 


(Accro ct al. 2010) 


PKS 2155-304* 


HBL 


0.117 


15.7 


63 


HESS 


3.34 ± 0.05 Bt at ± 0.1 syst 


14% 


400 


(Aharonian ct al. 2009) 












low 


3.53 ± 0.06 stat ± 0.10 syst 


15% 


200 


(Abramowski et al. 2010b) 


BL Laccrtac 


LBL 


0.069 


14.28 


35 


high 


3.64 ± 0.54 Btat ± 0.2 syst 


3% 


200 


(Albert et al. 2007c) 


1ES 2344+514* 


HBL 


0.044 


16.4 


10 


average 


2.95 ± 0.12 stat ± 0.2 syst 


11% 


200 


(Albert et al. 2007a) 


H 2356-309 


HBL 


0.165 


17.24 


8 


average 


3.06 ± 0.15 atat ± 0.10 syst 


2% 


240 


(Abramowski ct al. 2010a) 



a (Miller ct al. 1978; Lanzctta ct al. 1993) 



b IIFrandini et al.ll2(Till) 

c iMazin & Goebe]||2007l) 

* Blazars that are reported as variable in the TeV band, according to TeVCat i jhttp: //tevcat.uchicago.edu/} . 

TABLE 1 

GeV-TeV properties of the VHE blazar sample taken from the literature. Columns (1), (2) and (3) show the spectral 

ENERGY DISTRIBUTION (SED) TYPE, REDSHIFT , AND SYNCHROTRON PEAK FREQUENCY (v S yn), RESPECTIVELY. Fermi VARIABILITY INDICES 

(4) WERE TAKEN FROM THE 1FGL CATALOG (IAbDQ ET A L. 2010b). Fermi STATES (5) ARE IDENTIFIED IN THIS WORK USING 27 MONTH 
Fermi LIGHT CURVES AS DESCRIBED IN SECTION [5] In CASES WHERE Fermi DATA ARE CONTEMPORANEOUS WITH TeV OBSERVATIONS, THE 
CORRESPONDING TeV INSTRUMENTS ARE LISTED IN COLUMN (5). TeV SPECTRAL INDICES (6) WERE TAKEN FROM THE REFERENCES LISTED 
(9). TeV INTEGRAL FLUXES (7) ARE ABOVE THE LISTED ENERGY THRESHOLD (8) AND IN UNITS OF CRAB NE BULA FLUX. FOR THE CRAB 

Nebula unit conversions, spectral measurements above 350 GeV from (ICelikII2008I) are used. 



Twelve out of 26 blazars did not have enough statis- 
tics for a temporal separation of the Fermi data set into 
different flux states. Therefore for this subsample, an av- 
erage spectrum was calculated using the entire data set. 
Data from another subsample with 12 blazars were split 
into high and low-flux states as described in Section 
Data from the two brightest blazars (Markarian 421 and 
3C 279) were split into three subsets, with low, medium 
and high-flux states. See Table [2] for a summary of our 



Fermi data analysis results. 

O ur analysis results are consistent with the 2FGL cata- 
log (|Nolan et al.ll2012t) . We used the combined GeV- TeV 
SEDs (see Figure [8|) to estimate the IC peak frequency 
band of each blazar (see Section l6~Tj) . Our sample con- 
tains a handful of candidate "TeV-peaked" blazars that 
we discuss in Section 16.21 In addition, considering the 
fact that Fermi spectral indices do not vary significantly 
between low- and high-states, we studied the change in 
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Fig. 6. — TeV photon index distributions for our sample (see Ta- 
bic [TJ. For blazars that have multiple published results in Table [T] 
the most recent one was used. Top, middle and bottom panels 
show HBLs, IBLs+LBLs. and FSRQs, respectively. HBLs tend to 
have harder spectra than the rest of the sample. FSRQs have the 
softest spectra. Note that the TeV indices are not EBL-corrected. 

spectral index from GeV to TeV as a function of the red- 
shift, thus confirming the EBL effect on TeV spectra with 
a model- independent approach (see Section [6"T3l) . On the 
other hand, interesting spectral features in the GeV band 
are observed. To probe these features, the data were fit- 
ted with three different functions and the corresponding 
fit improvements were calculated (see Section |6T4")) . Fi- 
nally, in Section 16.51 we extended this study to contem- 
poraneous combined SEDs. 

6.1. IC Peak Frequency 

The peak frequency of the IC component is a salient 
parameter for describing blazar non-thermal continua 
and studying population trends. Systematic studies for 
measuring the IC peak frequency mostly suffer from the 
lack of statistics and simultaneous data . A similar work 
was carried out in (jZhang et al.l 120121 ) , where archival 
multiwavclcngth data were used to model TeV blazar 
SEDs and determine the IC peak frequency (i^c)- A 
positive correlation between v syn and vie was reported. 
In this work, we focus on finding the IC "peak frequency 
band" rather than the "peak frequency" , using a model 
independent approach. For each blazar SED shown in 
Figure [SJ we identify the energy decade in which the 
largest amount of power is emitted. Note that the spec- 
tral points used in the VHE spectra are EBL-corrected. 
Figure [9] shows the distribution of the IC peak bands for 



different blazar types. We observe that the FSRQs, LBLs 
and IBLs have the maximum of their emission mostly 
below 1 GeV. On the other hand, HBLs tend to peak 
in the TeV range. This positive correlation between the 
synchrotron (J4yn) an d the IC peak frequencies {v\c)i 
is in accordance with simple SSC models that predict 
a posi tive correlation between v syn and v\q ()Abdo et al.l 
l2010al) . The dashed lines represent the same distribu- 
tions with the bright AGN sample from t he first three 
months of Fermi data (Abdo et al. 2010a|). Our results 
tend to span the high frequency sides of all distributions 
and one clearly sees a shift to higher frequencies in the 
case of HBLs. This is expected since our sample consists 
of TeV-selected objects, that mostly correspond to rela- 
tively weak sources in the GeV data, and are therefore 
less likely to appear in a bright AGN sample. It should 
also be noted that we use a model indepen dent method 
using only Fermi and VHE data, whereas (jAbdo et al.l 
2010a) uses multiwavelength data and some modeling in 
cases where the soft X-ray band is dominated by the 
synchrotron component, a typical feature for our blazar 
sample. 

6.2. Hard TeV BL Lac Objects 

The combined GeV- TeV spectra of some blazars in 
our sample (1ES 0229+200, 1ES 0347-121, 1ES 1101- 
232, 1ES 1218+304, H 1426+428) suggest a v lc beyond 
~ 1 TeV. These blazars are mostly weak or non-detected 
in the Fermi range, with a hard spectral index in both 
GeV and TeV bands. It follows that they may belong 
to the so-called ultra- hig h- frequency-peaked BL Lac sub 
class (UHBLs; see, e.g., (jCostamante 11201 Tl) ) that would 
constitute the extreme end of the population, and is ex- 
pected to dominate the TeV luminosity of the universe. 
Several mechanisms have been set for th to explain th e 
formation of these hard 7-ray spectra (jLefa et al.H2011l) . 
Extensive spectral analysis of these objects would be 
valuable for EBL and intergalactic magnetic field mea- 
surements. It should be noted that at energies of a 
few TeV and beyond, our spectra become EBL-model- 
dependent. For this reason, we have compared our 
adopted EBL model with two other models from r ecent 
studies (jFinke et al.ll2010t iFranceschini et al.ll2008D . We 
have found that for the data samples mentioned above, if 
we used any of the other two EBL models, the dispersion 
in highest energy flux points would be less than 20%, and 
consequently the observed spectral upturns would not be 
affected significantly. 

a deeper variability study 
would relate to arguments 
ray production as the ori- 
gin of TeV blazar emission, since in that scenario no 
short ti mescale variability w ould be expected to be ob- 
served (jMurase et al.l 1201 21 ). Among the UHBL can- 
didates, the ones that are present in the 1FGL cata- 
log (1ES 1101-232, 1ES 1218+304, H 1426+428) have 
relatively small Fermi variability indices (see Table [lj . 
In addition to that, our calculations of F val for all five 
blazars using 27 months of Fermi data do not indicate a 
significant hint of variability either (see Table [5]) . 

6.3. Spectral Variability 

VHE emission from blazars is highly variable. This 
variability, manifested in irregular flares, is one of the 



With additional data 
carried on these blazars 
that support the cosmic 
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a jMiller et aUU978tlLanzetta et al.lll993H 
b l lPrandini et alj|2011l) 
c UMazin & Goebelll2007l) 

Blazars that are reported as variable in the TeV band, according to TeVCat( http://tcvcat.uchicago.edu/}. 



TABLE 2 

Fermi ANALYSIS RESULTS FOR THE SAMPLE. COLUMNS (1) AND (2) ARE THE SAME AS IN TABLE [T] F va , r (3) IS THE 
CALCULATED FLUX VARIABILITY AMPLITUDE (SEE SECTION 16. 3p FOR THE 27 MONTH PERIOD. Fermi STATES (4) ARE AS 
DESCRIBED IN SECTION In CASES WHERE Fermi DATA ARE CONTEMPORANEOUS WITH TeV OBSERVATIONS, THE 
CORRESPONDING TeV INSTRUMENTS ARE LISTED IN COLUMN (4). TggV (5) REPRESENTS THE PHOTON INDEX AND fi-ioo (6) 
THE INTEGRAL FLUX FOR 1-100 GeV. TEST STATISTICS (TS) AND LIVE TIME CORRESPONDING TO THE LISTED FLUX STATE 

ARE GIVEN IN COLUMNS (7) AND (8), RESPECTIVELY. 



most typical and promising blazar behaviors for study- 
ing the nature of underlying emission mechanisms. The 
observed flux chan ge during a VHE flare can be as rapid 
as minute scales (|Mkn 501: Albert et al.l l2007ef) and as 
large as 40 times the baseline emission ()Arlen et al.l 
I2013T) . Blazars that have been reported to have a vari- 
able flux are marked with an asterisk in Table [T] On the 
other hand, Fermi data do not exhibit flux variability as 
extreme as in the VHE band. In fact, having a smaller 
effective area than the ground-based VHE telescopes and 
operating mostly in survey mode rather than pointing, 
Fermi-LAT does not have the sensitivity to probe sub- 
hour timescale variability in blazars. Still, a possible 
corre lation between GeV and TeV emission remains vi- 
able (|Abdo et al.H2011b lAleksic et aLll2011b| ) and an en- 
hanced activity in the high-energy tail of the Fermi band 
could therefore indicate a TeV flare. In this frame, mon- 
itorin g GeV flares to trigge r TeV observations is impor- 
tant (jErrando et al.l [201 laf ). and potentially could help 
in probing fast variability. To examine variability within 
the Fermi data, we compared high- and low-state Fermi 
spectra from 14 blazars (see Table [3]). Half of these 
blazars have their integral flux in 1-100 GeV (-Fi_ioo) 
increased by at least 90% in the high state. The largest 



flux increase is seen in the case of the two FSRQs 3C 279 
and PKS 1222+21. As depicted by their respective light 
curves in Figure!?) these two objects have undergone dra- 
matic GeV flares. Such a large scale flux increase does 
not hold for the remainder of the blazars. However, one 
should keep in mind that for most of the TeV blazars, 
the Fermi band is a relatively stable region of the SED, 
since it samples the low energy part of the parent elec- 
trons, that have a longer cooling time. Table |3] gives a 
summary of the results of the spectral variations seen in 
the Fermi data. We also calculated the variability am- 
plitude (F var ) within 27 months of Ferm i data for each 
blaza r, using the method described in (jVaughan et al.l 
2003). F var is a measure of the intrinsic source vari- 
ance, calculated based on excess variance. For blazars 
with negative excess variance, 95% confidence level up- 
per limits are given. The blazars 3C 66A, PKS 1222+21 
and 3C 279 are the most variable ones according to this 
calculation (F v ^ r > 0.5). Our results arc in agreement 
with the 1FGL catalog (see Table [I]) . Comparing these 
results with the TeV variability flags, we do not find any 
obvious relation between GeV and TeV variabilities (see 
Table [2]). 

Within the Fermi energy range, blazars in our sam- 
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SED type 


Increase in _Fi_ioo (%) 
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(1) 


(2) 
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3C66A 


IBL 


95 


2 00 + 02 


1 93 + 02 


S5 0716+714 


LBL 


100 


2 20 + 03 


2 1 3 + 04 


1ES 1011+496 


HBL 


30 


1 90 ± 04 


1 97 ± 04 


Mrk 421 


HBL 


90 


1 86 ± 03 


1 79 ± 02 


1ES 121 8+304 


HBL 


60 


1 68 + 08 


1 73 + 1 


W Comae 


IBL 


110 


2 10 ± 06 

— ■ 1 1/ — 1 — UiUU 


2 07 ± 06 


PKS 1222+21 


FSRQ 


1290 


2 50 ± 04 

— ■ ' r V i — i — v / . v / l 


2.32 ± 0.02 


3C 279 


FSRO 


395 


2 54 + 03 


2 37 + 02 


PKS 1424+240 


IBL 


45 


1.82 ± 0.04 


1.87 ± 0.03 


PC 1553+113 


HBL 


30 


1.70 + 0.03 


1.74 + 0.03 


Mrk 501 


HBL 


80 


1.84 + 0.05 


1.83 + 0.04 


1ES 1959+650 


HBL 


60 


2.04 ± 0.06 


2.07 + 0.05 


PKS 2155-304 


HBL 


80 


1.95 + 0.03 


1.92 + 0.02 


BL Lacertae 


LBL 


140 


2.46 ± 0.05 


2.34 ± 0.04 



TABLE 3 

Spectral variations in the Fermi data for blazars for which at least two different Fermi flux states are 

AVAILABLE. COLUMN (1) SHOWS THE SPECTRAL ENERGY DISTRIBUTION (SED) TYPE. Fi_ioo IS THE INTEGRAL FLUX FOR THE 

1-100 GeV band. Column (2) shows the percent increase in Fi_ioo from low to high Fermi state. FSRQs and 
LBLs seem to show the most significant flux variability in this energy range. Columns (3) and (4) list the 
GeV photon indices in low and high Fermi states, respectively. No significant changes in photon index are 

SEEN, EXCEPT FOR THE TWO FSRQS FOR WHICH THE INDEX SHOWS A SLIGHT HARDENING FROM LOW TO HIGH FLUX STATES. 



pie do not exhibit dramatic changes in their spectral in- 
dex between different flux states (see Table [3|) . Conse- 
quently, this makes the photon index a reasonable pa- 
rameter to use for studying the non-contemporaneous 
combined SEDs. Figure [TU] shows a scatter plot of ob- 
served TtoV — TgoV versus redshift. A constant func- 
tion does not provide a good description for the data, 
with x 2 /dof = 204/27, which could be interpreted as 
a model-independent indication for the EBL absorp- 
tion. The difference between TeV and GeV photon in- 
dices increases with redshift. This is expected since the 
VHE 7-ray photons pair pro duce with the EBL pho- 
tons ( Fra nceschini et ail 120081 ) and this effect becomes 
more enhanced at larger redshifts, making the universe 
opaque to TeV 7 rays at distances larger than z ~ 0.5. 
HE spectra are not affected by the EBL, whereas VHE 
spectra become softer with increa sing redshift. A sim- 
ilar observation was reported by (|Abdo et al.ll2009t ). in 
a study carried out on a sample of TeV-selected AGNs 
detected with Fermi. 

Figure QT] shows the relation between the spectral in- 
dex TgcV an d the flux normalization _Fi_ioo obtained 
from power-law fits. FSRQs and two subgroups of BL 
Lacs are clearly separated in the parameter space. This 
is in accordance with the aforementioned positive corre- 
lation trend between v syn and v\c 1 since 1 GeV typically 
corresponds to the rising edge of the IC component in 
an HBL SED, sampling a relatively low flux with hard 
spectral index. On the other hand for an FSRQ, 1 GeV 
will correspond to the peak or the falling edge of the IC 
component. The fact that FSRQs have relatively more 
luminous IC emission explains the softening trend with 
a larger normalization factor. However, the pattern that 
we observe between different flux states of a given blazar 
is the opposite. In most slight spectral hardening 

accompanies high flux states, indicating a change in the 
spectral shape and enhanced flux increase at high-energy 
tail of the spectrum. 



6.4. Spectral Features 

In most of the blazars in our sample, we observe inter- 
esting spectral features in the Fermi band, that appear 
as dips in the 1-100 GeV energy range. In an attempt 
to find a quantitative description for these features, we 
fit the Fermi spectral points with a simple power law 
(PL; Equation @) and a broken power law (BPL; Equa- 
tion and then compare the results. 



dN/dE = N x 



f(E/E b ) 
\(E/E b ) 



-r 2 



E < E b 
otherwise 



(3) 



In the PL fit, the normalization Nq and the spectral 
index T are free parameters, and the energy Eq is fixed 
at 1 GeV. In the BPL fit, the break energy E^ and the 
indices Ti and T2, along with the normalization Nq are 
free. In Table|3J we list the best-fit parameters from both 
functions and the likelihood ratio test results of BPL over 
PL. In 9 out of 33 cases, BPL yields a better fit over PL 
with more than 2a significance. 

There are several possible mechanisms that may cause 
the observed features in the SEDs. One possibility is a 
break in the electron spectrum caused by the synchrotron 
cooling effec ts, generally yielding a ch ange in spectral in- 
dex by 0.5 (jBottcher fc Chiandl200l) . which is in agree- 
ment with our results (see Table EJ. Another mecha- 
nism that could explain the observe d breaks is the ab- 
sorpt ion by an external photon field (jPoutanen fc Sternl 
2010). For those nine data sets where the BPL gives 
a better fit than the PL, the break energy ranges from 
~ 2 GeV to ~ 8 GeV. In addition, 7 of these data sets be- 
long to non-HBL blazars, that are usually characterized 
by broad emission lines, thought to be originating from a 
region of molecular gas (broad line region; BLR) that is 
highly ionized by the optically thin accretion disk. This 
seems in accordance with the idea of relating the Fermi 
spectral features to absorption of GeV photons on radia- 
tion from Hi (13.6 eV) and Hen (54.4 eV) recombination 
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Fig. 7. — Fermi-LAT aperture photometry light curves with no background subtraction, normalized to arbitrary units. The solid lines, 
representing the average counts per area per time, separate "low" and "high" flux states, that are later on used to produce "low" and 
"high" state spectra. In the case of Markarian421 and 3C 279 light curves, the dashed lines represent 1 a deviation from the average, 
dividing the data set into three separate flux states ("low", "medium", and "high"). The shaded areas show the contemporaneous time 
windows with the corresponding TeV instruments. 



continua in the BLR, that are expected to cause jumps 
in 7-r ay opacity around ~ 19.2 and ~ 4.8 GeV, respec- 
tively (jPoutanen k Sternll2010ft . We tested a general ab- 
sorbed power-law (APL) function of the following form 
on the Fermi data: 



dN/dE = No(E/E )~ 



T (_E,z,i5 abs ) 



(4) 



where the free parameters are the normalization N 
at Eq = lGcV, photon index T, and absorption line 
energy E a i, s . t 17 is the optical depth for the 7-7 pair 
annihilation of photons with energies E and -E a bs at a 
rcdshift of z. Within the Fermi energy band, BPL and 
APL functions fit the data equally well. Upturns at high- 
energy tails of Fermi spectra are observed (see, e.g., W 
Comae in Figure [5]) , but they are not statistically signif- 
icant enough to favor an absorption scenario over a BPL 
fit. Therefore, it appears that one cannot statistically 
distinguish between the BPL and APL fits, but possi- 
ble absorption scenarios are worth investigating further. 
To address this issue, we make use of contemporaneous 
GeV- TeV spectra to test and compare BPL and four dif- 
ferent APL scenarios (see Section l6~5j) . This permits us 
to test the APL over a larger energy range and investi- 
gate the apparent Fermi spectral absorption-like features 
with higher statistics. 

Another caveat related to these spectral features is that 
the upturn seen at the highest Fermi energy bin might 
be coming from a group of photons clustered in time. In 
that case the dip would be an artifact of a flaring event, 
thus not representative of the time-averaged spectrum. 
To make sure this is not the case, we checked the arrival 
times of the highest-energy photons and did not find any 
obvious clustering (see Figure [12")) . Note that the arrival 
time distributions should be considered within a given 
flux state. For instance, in the left panel of the figure, 
the red triangles represent the high energy photons from 
the high flux state and are evenly distributed in a time 
window that belongs to the high state. Therefore, one 
concludes that no clustering is found. 

6.5. Quasi- simultaneous GeV -TeV spectra 

Seven of the TeV spectra in our sample are contem- 
poraneous with Fermi observations and therefore merit 



a deeper analysis. We extended the work described in 
section T6.4I to this subsample. This time, in addition to 
PL and BPL fits, we tested four different scenarios of ab- 
sorption due to photons emitted from the BLR: Hi line 
(13.6 eV), Hen line (54.4 eV), H T k Hen combined, an d 
full BLR spectrum taken from (jPoutanen k Sternfe OlO). 

For single- and double-line absorption scenarios, Hi 
and Hen recombination continua are the most plausi- 
ble cases given that they are the most dominant ones 
in the BLR spectrum, and that the breaks we see in 
the Fermi spectra are located around a few GeV. As for 
the f ull BLR spectrum, taken from (jPoutanen k Sternl 
2010), it was modeled assuming a photoionized gas with 
the ionization parameter and the cloud density chang- 
ing with the distance to t he central ionizing source. 
See (jPoutanen k Sternll2010( ) for a detailed discussion on 
the 7-ray absorption within the BLR in the Fermi spec- 
tra. No general trend can be seen in the contemporane- 
ous data sample. BPL and full BLR absorption scenarios 
seem to fit well the combined spectra of the blazars PKS 
1424+240 and PKS 2155-304. A BPL (full BLR absorp- 
tion) function is preferred over the PL for the blazars 
PKS 1424+240 and PKS 2155-304 with a significance of 
~ 5<r (~ 4.8<r) and ~ 12er (~ 8.5<r), respectively (see Fig- 
ure EH]). The x 2 /dof values of PL, BPL and APL fits are 
32/9, 3.5/7, 5.3/7 for PKS 1424+240 and 148/8, 5.8/6, 
and 71/6 for PKS 2155-304, respectively. BPL fits yield 
Ar = 1.4 (PKS 1424+240) and Ar = 0.7 (PKS 2155- 
304), both larger than what electron cooling would pre- 
dict, which might indicate that an additional mechanism 
is at work. Both BPL and full BLR absorption scenarios 
provide a slight improvement in the MAGIC and VERI- 
TAS spectra of 3C 66A, albeit not significant. Similarly, 
for PKS 1222+21, BPL, Hi single line and Hi + Hen dou- 
ble line absorptions slightly improve the fit over PL. In 
the case of RGB J0710+091 and 1ES 1218+304, we don't 
observe any preference over the power-law fit. In case a 
7-7 absorption from BLR is at work, the cascades initi- 
ated in th is process might produce observ able GeV 7-ray 
emission (jRoustazadeh k Bo ttchcr 20TTI). and their syn- 
chrotron emission coul d contribute to the big blue bum p 
seen in several blazars (Roustazadch & Bdttchcr 2012). 
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Fig. 8. — GeV-TeV spectra for the sample of blazars in this study. Filled circles represent the Fermi spectra and the filled squares the 
TeV spectra. Considering the TeV flux state information given in the TeV papers, the best matching GeV and TeV spectral points are 
used for the combined analysis (shown in black). When available, spectral points belonging to other flux states (in both bands) arc plotted 
in gray. 3C66A (VERITAS and MAGIC respectively), RGB J0710+591, 1ES 1218+304, PKS 1222+21 (4C +21.35), PKS 1424+240, and 
PKS 2155-304 spectra arc quasi-simultaneous. 
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iV(xlO~ n ) 
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Figcv(x10- 12 ) 


ri 


r 2 


£ brcak (GeV) 


C"BPL 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


3C66A (low) 


2.78 + 0.07 


2.01 ± 0.02 


9.30 + 4.94 


1.81 + 0.06 


2.22 ± 0.07 


1.92 + 0.51 


4.12 


S5 0716+714 (low) 


1.85 + 0.06 


2.23 ± 0.03 


5.18 + 2.07 


2.14 + 0.04 


2.74 + 0.23 


5.43 ± 1.01 


3.07 


S5 0716+714 (high) 


3.68 + 0.13 


2.10 + 0.03 


2.29 ± 1.78 


2.00 + 0.06 


2.57 + 0.30 


4.10 ± 1.47 


2.29 


1ES 1011+496 (low) 


0.88 + 0.04 


1.94 ± 0.04 


1.70 ± 1.34 


1.70 + 0.10 


2.28 ± 0.20 


2.76 ± 1.15 


2.76 


PKS 1222+21 (low) 


1.72 + 0.06 


2.53 ± 0.04 


0.99 ± 0.70 


2.41 + 0.06 


3.11 ± 0.28 


3.36 ± 0.87 


2.84 


3C 279 (low) 


2.55 + 0.08 


2.59 ± 0.03 


0.15 + 0.28 


2.53 + 0.04 


3.72 ± 1.54 


7.70 ± 5.57 


2.51 


3C 279 (high) 


1.16 + 0.02 


2.39 ± 0.02 


8.14 + 10.72 


2.31 + 0.03 


2.74 + 0.26 


3.21 ± 1.79 


3.11 


PKS 2155-304 (low) 


4.19 + 0.10 


1.92 ± 0.02 


1.75 + 1.98 


1.86 + 0.03 


2.13 ± 0.14 


5.54 + 3.23 


2.08 


BL Lacertae (high) 


3.26 + 0.11 


2.39 ± 0.04 


1.08 + 0.09 


2.30 + 0.04 


3.00 ± 0.30 


4.49 ± 0.03 


2.50 



TABLE 4 

Fit results for power law (PL) and broken power law (BPL), where BPL yields a better fit over PL with more than 2a significance (9 out of 
33 cases). Columns (1) and (2) show the PL parameters, flux normalization at 1 GeV and the photon index respectively. Figcv (3) is the flux 

NORMALIZATION AT 1 GeV FOR BPL. N AND -FlGoV ARE IN ERG CM -2 S -1 . COLUMNS (4) AND (5) SHOW THE PHOTON INDICES FOR BPL, AS GIVEN IN EQUATION ([3J. THE 

break energy for bpl is listed in column (6). (xbpl (7) is the likelihood ratio test results of bpl over pl. for these 9 cases, the break energy -ebreak ranges 

from ~ 2 GeV to ~ 8 GeV. In addition, 7 of these data sets belong to non-HBL blazars. 
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Fig. 9. — Distribution of the IC peak bands, defined as the en- 
ergy decade in which the largest amount of power is emitted, in 
combined GeV-TeV spectra. Top, middle and bottom panels show 
HBLs, IBLs+LBLs, and FSRQs respectively. HBLs tend to peak at 
higher frequencies, in accordance with their respective synchrotron 
peak frequencies i> Byn (see Table [TJ, and a decreasing trend in 
IC peak bands from top to bottom panel is seen. The dashed 
lines represent the same distribution for the blazar sample from 
HAbdo et al.N2010al'l . 
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Fig. 11. — TgoV vs - nux m the energy band 1-100 GeV from 
the analysis of 27 month Fermi data with a power law fit (see 
Table El. Empty crosses, triangles and inverse triangles represent 
HBLs, IBLs+LBLs, and FSRQs, respectively. Solid lines connect 
different states of the same blazar. 




Fig. 10. — Ar (rxov — Tgov) vs - redshift. Empty crosses, 
triangles and inverse triangles represent HBLs, IBLs+LBLs, and 
FSRQs, respectively. The fact that a constant function does not 
provide a good description for the data could be interpreted as a 
model independent indication for the EBL absorption. 
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Fig. 12. — Arrival times of the photons in the highest energy bin (top panels) and aperture light curves with arbitrary flux units (bottom 
panels) for the blazars 3C 279 (left) and PKS 2155-304 (right). Blue upside-down triangles represent the low-state photons and red triangles 
the high-state ones. In both cases, the highest energy photons do not show any obvious clustering within their respective data sets. As 
described in Section [5] low- and high-states are distinguished based on the flux averages (solid lines) in light curves. 
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Fig. 13. — Contemporaneous GeV-TeV spectra with power-law (dashed lines), broken power-law (solid lines) and power-law with full- 
BLR-absorption fits. BPL and full BLR absorption scenarios seem to fit well the combined spectra of the blazars PKS 1424+240 and 
PKS 2155-304. A BPL (full BLR absorption) function is preferred over the PL for the blazars PKS 1424+240 and PKS 2155-304 with a 
significance of ~ 5<r (~ 4.8a - ) and ~ 12(X (~ 8.5c) respectively. 
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7. SUMMARY 

We study blazar spectral properties with a focus on 
the GeV-TeV energy range for a sample of VHE blazars. 
In order to obtain a set of joint GeV-TeV blazar spec- 
tra, we analyze the first 27 month Fermi data for VHE 
blazars and combine our results with archival VHE data. 
In cases where the Fermi data set does not overlap with 
the TeV observations but has enough statistics, we split 
the data into high and low flux states and join the best- 
matching subset with the corresponding TeV spectrum. 
The peak frequency band of the inverse Compton com- 
ponent increases following the order FSRQ -> LBL&IBL 
-> HBL. Thus, our results confirm the positive correla- 
tion between ^ syn and z^ic- We note that Fermi spec- 
tra from different flux states for a given TeV blazar do 
not undergo a significant change in photon index. The 
variability amplitudes within our Fermi data set do not 
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show an immediate correlation with the reported TeV 
variabilities for individual blazars. We find that in many 
cases a power law is insufficient to describe the GeV-TeV 
spectra and a broken power law improves the fits, espe- 
cially for non-HBL blazars, where the BLR emission may 
have an effect on the observed spectral shape. In some 
blazars we observe absorption-like spectral features. We 
present seven quasi-simultaneous joint spectra, on which 
we test possible absorption scenarios from the BLR. Even 
though the absorption seems to describe well some of the 
observed spectra, no general pattern can be identified. 
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